Matter-wave localization in a weakly perturbed optical lattice 
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By numerical solution and variational approximation of the Gross-Pit aevskii equation, we studied 
the localization of a noninteracting and weakly-interacting Bose-Einstein condensate in a weakly 
perturbed optical lattice in one and three dimensions. The perturbation achieved through a weak 
delocalizing expulsive or a linear potential as well as a weak localizing harmonic potential removes 
the periodicity of the optical lattice and leads to localization. We also studied some dynamics of 
the localized state confirming its stability. 
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I. INTRODUCTION 

Anderson suggested the possibility of localization of an 
electron wave in a nonperiodic random potential about 
fifty years ago pQ. After the recent experimental real- 
ization of Anderson localization of a Bose-Einstein con- 
densate (BEC) [2j [3], this topic has been the subject of 
intense research activities. Billy et al [2 observed the 
exponential tail of the spatial density distribution of a 
87 Rb BEC after releasing it into a one-dimensional (ID) 
waveguide with a controlled disorder potential created 
by a laser speckle. Roati et al. [3] observed the Ander- 
son localization of a noninteracting 39 K BEC in a ID 
quasiperiodic bichromatic optical lattice (OL). Recently, 
the experimental realizations of three-dimensional (3D) 
localization of ultracold gases of 40 K [4] and 87 Rb [5] 
atoms in a 3D speckle potential were also reported. In 
the experimental studies of localization of a BEC, bichro- 
matic OL [3] [6] , shaken OL [7 , cold atom lattice [8] and 
speckle potentials [2] have been used. There have also 
been studies of the effect of a repulsive atomic interac- 
tion on localization [9]. Bichromatic OL [T0HT2] , cold 
atom lattice [13] [14] , random potentials [T5UTF , and OL 
plus harmonic trap [19], among others, have been used in 
theoretical studies. Quasiperiodic [I9j[20] or random [2TT - 
[24] potentials for Anderson localization were realized by 
optical means. In Ref. [T9] . localization and Bloch oscil- 
lation of a BEC was studied in two spatial dimensions for 
OL perturbed by a harmonic potential [25] . Here we shall 
consider localization in ID and 3D in the same potential. 

However, there are other ways of obtaining a nonpe- 
riodic OL. Here we consider another approach for per- 
turbing the OL, so that the periodicity of the original 
OL is completely removed. If we perturb the OL by a 
weakly delocalizing (i) linear potential or (ii) expulsive 
harmonic potential, the periodicity of the OL is removed 
without adding an extra confining term to the OL. Both 
linear [26j [27] and expulsive [28] harmonic potentials can 
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be made in laboratory by optical and magnetic means. 
In Ref. [27] a tilted trap was used for the observation 
of Bloch oscillations in a repulsive BEC, where the ex- 
ternal potential is simply given by the sum of the linear 
(gravitational) component and the OL laser field. In Ref. 
[28] an expulsive potential was used to study solitons in 
BEC. 

Here, with the numerical solution of the Gross- 
Pitaevskii (GP) equation, we investigate the localization 
of cigar-shaped BECs in a ID OL, weakly perturbed by 
linear and expulsive confining harmonic potentials. We 
also considered adding a weak confining harmonic trap 
to the OL, which removes the periodicity of the OL to 
obtain localization. The localized states so obtained are 
completely distinct and had a different spatial extension 
from the states trapped in the confining harmonic trap 
alone and hence were localized due to the nonperiodic 
OL. It is found that the localized states in the three types 
of periodicity-broken OL are practically the same and do 
not depend on the weak perturbation used to remove the 
periodicity of the OL. In most cases the density profile 
of the central part of the localized state is quite similar 
to a Gaussian shape. Hence it is worthwhile to solve the 
GP equation analytically using a Gaussian variational 
ansatz. Indeed we find that the variational solution fits 
well the density distribution in the central region, espe- 
cially when the localization takes place essentially on a 
single site of the perturbed OL. However, the localized 
states usually have a long undulating exponentially de- 
caying tail in either direction and extends to several sites 
of the OL corresponding to Anderson localization in a 
weak potential p~5j [24]. The stability of the localized 
state under small perturbation is also investigated. 

We also extend our study to localization in 3D. In 3D, 
the effect of disorder on an initial state at large times 
depends strongly on the energy, and a mobility edge 
was predicted, which corresponds to a critical energy at 
which a transition from localized to extended character 
of the state occurs [2^ . As observed in the recent experi- 
ments in 3D [4, 5 , a two-component density distribution 
emerges at large times consisting of an expanding mobile 
component and a nondiffusing localized component after 
the ultracold gas of atoms are released into an optical 
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speckle field. In this paper, different from the experi- 
mental scenario, we focus our attention on the stationary 
localized states, and not on the large-time expansion dy- 
namics. We find that the weakly perturbed 3D OL also 
leads to localization. In the 3D realm, we only present 
results for the expulsive perturbation of the OL, as other 
perturbations lead to similar results. In this case we show 
results for noninteracting BEC as well as for weakly re- 
pulsive BEC. The numerical chemical potentials are in 
good agreement with the variational results in all cases. 
The localized states of the 3D BEC are found to be stable 
against small perturbations. 

In Sec. [IT] we present the ID and 3D models as well 
as their variational analysis used in this study of local- 
ization. In Sec. |III| we study numerically the statics and 
dynamics of a ID and a 3D BEC by analyzing the atom 
density profiles and chemical potential of the localized 
states in a perturbed OL, and compare the numerical 
results with the variational solutions. The study of a 
sustained small oscillation of the localized states upon 
external perturbation confirmed their stability. In Sec. 
IV we present a summary of our findings. 



II. ANALYTICAL CONSIDERATION 

The dynamic behavior of a weakly interacting Bose gas 
at low temperature on a 3D OL is well described by the 
GP equation which has the form of a (3+l)-dimensional 
nonlinear Schrodinger equation [30] : 

*W = -\v 2 $ + G\$\ 2 $ + U(x,y : z)$, (1) 
U(x,y,z) = -Vo[cos(2x) + cos(2?/) + cos(2z)], (2) 

where = $(x,y,z,t) is the BEC wave function with 
normalization / // dxdydz\&\ 2 = 1, and V 2 is the 3D 
Laplacian. The constant G = 4iraN characterizes the 
two-body interaction, TV the total number of atoms, a the 
atomic scattering length, and U the 3D OL. The spatial 
variable, time, density \&(x, y, z, t)| 2 , and potential are 
measured in units of Iq = \/(2ir), uiIq/H, Iq 3 and 2Er, 
respectively, where Ef> = h 2 / (2m\ 2 ) is the recoil energy 
of an atom absorbing a lattice photon, A is the wavelength 
of the OL potential, and m is the atomic mass. 

We consider a 3D trapping potential U (x, y, z) made 
up of three weakly perturbed OL components defined as: 



U(x,y,z) = ^2v(xi), 



(3) 



2=1 



where #1,2,3 correspond to x,y,z, respectively. Each of 
the three components has the following forms: 



V(xi) = -Vb cos(2xi) - Cxi, 
V(xi) = -V cos(2xi) ±Cx 2 , 



(4) 
(5) 



with Vb = 5, corresponding to perturbations by a weak 
linear [Eq. Q] and confining [Eq. ([5|, positive sign] and 
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FIG. 1: (Color online) The linear, confining, and expulsive 
harmonic potentials V(x) as a function of x in dimensionless 
units from Eqs. and Q for C = 0.005, V = 5. 



expulsive [Eq. ([5|, negative sign] traps, with strength C. 
These perturbations slightly modify the periodic OL and 
generate nonperiodic OL potentials. 

We also consider a cigar-shaped condensate oriented 
in the longitudinal direction (x direction). The trapping 
potential is a weakly perturbed OL in the longitudinal 
direction and the transverse (y, z directions) dynamics of 
the condensate is frozen to the respective ground states 
of harmonic traps. This cigar-shaped system will be de- 
scribed by using an effective ID model derived from the 
3D GP equation by integrating over the transverse 
variables [3T] , Then, the time-dependent wave function 
(j> = (/>(#, t) of a BEC can be described by the following 
ID GP equation [3U [32] 



dt 



ld 2 (/) 
'2dx~z 



•<#| 2 </> + F(*)</>, 



(6) 



with normalization \<j)\ 2 dx = 1, nonlinearity g = 
2aN/a 2 L , and potential V(x) given by Eq. Q or (|5]). 
The spatial variable x, time t, density |0| 2 , and energy 
are expressed in normalized units a± = y/H/ (mu±), cjJ 1 , 
a^ 1 and Huj±, where uj± is the angular frequency of the 
transverse trap. These weakly perturbed OLs along x 
axis are illustrated in Fig. [l] which shows that the OL 
dominates for small x, but the periodicity of the original 
OL is completely removed. 

Stationary solutions u(x,y,z) to Eq. ([T]) can be given 
by 2/, z, t) = exp(—ijjit)u(x : y, z), with /i the chemical 
potential. For convenience, we use unified notation u to 
denote the real wave function with u = u(x) for ID and 
u = u(x,y,z) for 3D BEC, respectively. We choose d(= 
1,3) to denote the dimension of space. The stationary 
real wave function obeys 



1 ^ f)2 d 



(7) 
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The nonlinearity G should be replaced by g in ID. Equa- 



3 



tion (|7| can be derived from the Lagrangian density 
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To apply the variational approximation to Eq. ([7]), we 
use the variational Gaussian ansatz: 



u = 



rV4 



^n-p(-i), (9) 

2=1 V 7 



with iu the width and Af the normalization of the local- 
ized BEC. The effective Lagrangian of the system can be 
found upon by substituting Eq. ([9| into Eq. (JsJ) and 
integrating over space [33J: 



L = fiAf — [i 



dM _ GM 2 
4uP 2~~ 
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dVoATe- 



- Rd 



Mw 2 



(10) 



where R = for the linear potential and R = ±C for the 
confining and expulsive harmonic potential. The Euler- 
Lagrange equation dL/dfi = yields the normalization 
J\f = 1. We use it in the following equations. The remain- 
ing equations dL/dw = dL/dJV = yield, respectively, 



1 



G 



W Q w A exp(-w 2 ) + 2#w 4 ,(il) 



G 



, o ■ r -dV exp(-w 2 ) + ,(12) 

4w 2 (wV2^) d % 



and determine the width w and the chemical potential 
fi. The linear perturbation has no effect on the width 
and chemical potential of the stationary localized states 
because R = for the linear perturbation. 



III. NUMERICAL RESULTS 

We performed the numerical integration of GP equa- 
tions and ([6| employing the imaginary or real-time 
split-step Fourier spectral method with space step 0.04, 
time step 0.001. We checked the accuracy of the results 
by varying the space and time steps and the total num- 
ber of space and time steps. We also checked the results 
with those obtained from a split-step Crank-Nicolson al- 
gorithm using same steps [34] . 

By numerically solving equation ([6]), we first investi- 
gate the effects of the weakly perturbed OL on a lo- 
calized state in ID BEC. To obtain a stationary local- 
ized state, the initial input pulse is taken as 0(x, 0) = 
exp(— x 2 /2)/7r 1 / 4 with a parabolic trap V'{x) = x 2 /2 in 
our numerical integration. During the numerical simula- 
tion, the parabolic trap is slowly turned off, and the po- 
tential Q or ([5| and the nonlinearity g is slowly turned 
on with a very small increment (= 0.00001) in each time 
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FIG. 2: (Color online) (a) Normalized numerical density pro- 
file |^(x)| 2 of the localized BEC versus normalized x for po- 
tential @ with C = 0.1,0.05, and 0.005 and g = 0. (b) 
Normalized numerical (Num) and variational (Var) chemical 
potentials \i of the localized BEC with g = versus C for the 
confining and expulsive perturbed OL 



step. The time evolution is continued until the parabolic 
trap is completely turned off and the potential Q or ([5| 
is completely turned on, and the the nonlinearity g has 
the desired value. 

In Figs. [2] we exhibit typical numerical and variational 
results of the localized states with different perturbations 
for g = 0. In the case of the linearly perturbed OL [Eq. 
Q], Fig. [2] (a) shows that the numerical density profiles 
of the central part of the localized states are the same for 
C = 0.005,0.05 and 0.1, and the numerical chemical po- 
tentials are all —2.897 in agreement with the variational 
result, —2.894. For confining and expulsive perturbed 
OL ([5]), in Fig. [2] (b) we plot the numerical and varia- 
tional chemical potentials \i of the localized states. In the 
case of the expulsive perturbation, the chemical poten- 
tial decreases with the perturbation strength. The oppo- 
site happens for the confining perturbation. The reason 
is that an expulsive perturbation weakens the trapping 
of the OL and a confining perturbation strengthens the 
trapping. For small C, the influence of C on the density 
profile and chemical potential of a localized state is small. 
We also find that the stability of the localized states is 
related to the perturbation strength C. If we increase C, 
the stability situation does not improve. But if we reduce 
C the localized states could be unstable as C — >• 0. So 
we keep C > 0.005 in this paper. 

With the periodic potential — 5cos(2x), the noninter- 
acting Schrodinger equation [viz. g = in Eq. ([6|] per- 
mits only delocalized states in the form of Bloch waves. 
Localization is possible in the noninteracting Schrodinger 
equation due to the nonperiodic nature of the perturbed 
OL, as shown in Figs. [3] (a) and (c) where we exhibit 
the typical density profiles of the localized states in the 
linear and confining harmonic potentials for g = and 10 
and compare them with the corresponding variational re- 
sults. To confirm that these localized states are induced 
by the weak perturbation of the OL, in Fig. [3] (c), we 
also show the localized state for g = for the confining 
harmonic potential alone (see the black dot-dashed line). 
This state has a different spatial extension and can be 
distinguished from the other states localized due to the 
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FIG. 3: (Color online) Numerical (Num) and variational (Var) 
density profile |^x(x)| 2 of the localized BEC versus x for the 
(a) linear and (c) confining harmonic perturbations of the OL 
with C = 0.005, and g = and 10, respectively. In (c), we 
also show the density profile of the localized BEC trapped by 
only the confining harmonic potential V(x) — 0.005x 2 with 
g = (the black dot-dashed line). The numerical (Num) and 
variational (Var) chemical potentials fi of the localized BEC 
versus g for (b) linear and (d) confining perturbations of OL. 
All quantities are dimensionless. 
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FIG. 4: (Color online) Center of the localized state xo ver- 
sus time t during the oscillation of the localized BEC on the 
nonperiodic OL Q and (J5| initiated suddenly by introducing 
an initial momentum po = 0.1 through the transformation 
u(x) — >> uo (x) exp(ip x) with g = and C = 0.005 for the (a) 
linear, (b) confining, and (c) expulsive perturbation, respec- 
tively. 



perturbation of the OL. Figures [3] (b) and (d) show, re- 
spectively, the effect of the nonlinearity on the chemical 
potential (again in the agreement with the prediction of 
the variational results). The corresponding results for 
the expulsive perturbed potential are in agreement with 
those of the confining potential, and are not shown here. 
To understand the effects of g on /i, and on localization in 
general, it is useful to consider the variational equations 
(11) and (12); for g = 0, fi = —2.894. As g is increased, 
H increases and the system is more repulsive and extends 
over more region in space and the localization length in- 
creases. Eventually, at g « 18.6 the state is completely 
delocalized and beyond this value of g, no localization is 
possible and there is no real solution of Eq. (11) for the 
width. 

As a typical Hamiltonian system, the stability of an or- 
dinary localized BEC can be determined by the Vakhitov- 
Kolokolov (VK) criterion [35]. Of course, the direct ap- 
plicability of this criterion depends on the form of the 
nonlinearity and on the form of trapping potential [36] . 
In the case of an attractive BEC confined in a mag- 
netic trap, for example, a general stability criterion is 
derived for the ground states of the GP equation in Ref. 
[37] . This criterion states that a localized state is sta- 
ble only if the chemical potential \i and the number 
of atoms N = \u(x)\ 2 dx satisfy dN/dfi < 0[37]. 
The chemical potential \i corresponds to —A of Ref. 
[37] . As for the repulsive nonlinearity, however, the sta- 
bility of gap soliton families obeys an inverted ("anti- 



VK") criterion dN/dfi > [38]. In Eq. the num- 
ber of atoms N is linearly related to the nonlinearity 
via g = 2aN/a 2 \_, and the wave function is normalized 
J_ oo \u(x)\ 2 dx = 1. Then, the anti-VK criterion can be 
expressed as dg/dfi > 0. Figures [3] (b) and 3(d) show 
that the anti-VK criterion dgjd[i > holds, hence the 
localized states are stable in the parameter range. To jus- 
tify the anti-VK criterion, the stability of the localized 
BECs is tested in systematic numerical integration of GP 
equation |6| with g > 0. For example, with g = 10, first 
we create a localized BEC in a nonperiodic OL. Then, 
the stability of the localized states was established by 
slightly changing the potential (by multiplying by a fac- 
tor of 0.98) in the real-time routine and continuing the 
time evolution to 10000 units of time. We find that the 
localized states remain localized in real-time routine and 
hence are considered to be stable. At large time scales, 
the behavior of an initially localized wave packet has been 
analyzed numerically in Ref. [18 . It was shown that 
above a certain critical strength of nonlinearity a subdif- 
fusion occurs. We emphasize here that, in our numerical 
calculations, the stationary localized states are obtained 
in an adiabatic way, and the stability of the localized 
states was tested within a short period of time. 

Nevertheless, the anti-VK criterion can not be natu- 
rally extended to the noninteracting BEC (viz. g = 0). 
To test the stability of these noninteracting localized 
BEC, first we create a localized BEC in a nonperiodic 
OL. Successively, at t = 0, we suddenly introduce a 
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FIG. 5: (Color online) (a) Normalized numerical density pro- 
file \u(x)\ 2 of the localized BEC versus x with g — and 
C = 0.005 for the linear, confining and expulsive perturba- 
tions, respectively, (b) Normalized numerical density pro- 
file \u(x)\ 2 of the localized BEC versus x with g — and 
C = 0.005, 0.05, 0.1 for the linear perturbations. 



small initial momentum po = 0.1 by changing u(x) — >• 
uq(x) ex.p(ip x), where uq(x) is the wave function of the 
localized BEC at t = 0. The center xq of the localized 
states versus time t is presented in Figs. [4]where xq is ob- 
tained by instant Gaussian-function fitting to the central 
region of the density distribution. It is found that after 
the perturbation, the localized states perform a sustained 
oscillations but remain localized which imply that the 
localized states are stable against small perturbations. 
Contrary to our intuition, however, the oscillations are 
quasiperiodic because of the deformation of the density 
profile, as explained in Ref. [39] . 

One interesting earmark of Anderson localization in a 
weak disordered potential is a long exponential tail of 
the localized state inspite of the Gaussian distribution of 
the central part [40]. To observe the effect of different 
perturbations of the OL on the tail of localized state, 
we present in Fig. [5] the density profiles of the localized 
BEC on log scale. As shown in Fig. [5] (a), the localized 
BECs in the three perturbed OLs have long exponential 
tails extending far beyond the central Gaussian distribu- 
tion. Further investigations show that the exponential 
tails are symmetric with respect to x = for the confin- 
ing and expulsive perturbations, i.e., both left and right 
localization lengths are the same. In these cases, the lo- 
calization length decreases with the increase of C. For 
linear perturbation, both the perturbed OL and the ex- 
ponential tails are asymmetric with respect to x = 0, and 
the asymmetry increases with the increase of the pertur- 
bation strength C, as shown in Fig. [5](b). Although, the 
localization lengths depend on the perturbation strength 
C, the dependence is weak for small C. 

Now we study localization in the 3D system with a 
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FIG. 6: (Color online) The 3D contour plot of normalized den- 
sity profile of the localized state for G — from (a) numerical 
calculation and (b) variational approximation with the 3D OL 
perturbed by the expulsive potential ([5} with C = 0.005. In 
the plot the contour was set at a density of 0.001. 
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FIG. 7: (Color online) (a) and (b) The same as Figs. 6 for 
G = 10. (c) Chemical potential fi for different G from nu- 
merical simulation and variational approximation, (d) Width 
w for G = 10 during real time evolution when at t — 25 the 
nonlinear ity G was suddenly changed from 10 to 11. 



perturbed OL Eq. (J3j) . As in a ID system, the results for 
three types of perturbations are very similar in the 3D 
system. Hence, in a 3D system we illustrate the results 
for the expulsive perturbation only. First, we illustrate 
the results for the noninteracting case: G = 0. To test 
the stability of the localized state for G = we use the 
real-time propagation of the 3D GP equation with 
the initial state at t = taken as the variational local- 
ized state. The initial solution settles down to the final 
localized state after iteration during an adequate interval 
of time. The 3D contour plot of the density profiles of 
the localized state with its initial variational counterpart 



6 



are shown in Figs. [6] (a) and (b). Although the varia- 
tional solution is spherically symmetric, this symmetry 
is absent in the actual localized state. This is because 
the 3D perturbed OL is not spherically symmetric. In 
spite of this, the variational approximation is very useful 
in the study of localization, and the numerical solution is 
nearly spherically symmetric. The variational chemical 
potential is —8.6826 compared to the numerical result of 
-8.690. 

Next, we study 3D localization in a repulsive system 
with G = 10. Numerically, we find the stable localized 
state by real-time propagation with the variational ap- 
proximation as the initial input. The 3D contour plots 
of the localized state with its variational counterpart are 
shown in Figs. [7] (a) and (b). For G = 10, the nu- 
merical profile of the localized state in Fig. [7] (a) shows 
more deviation from spherical symmetry compared to [6] 
(a) for G = 0. The variational chemical potential of 
the localized state in Fig. [7] (a) is —5.604 compared to 
the numerical result of —5.65. In Fig. [7] (c) we show 
the numerical chemical potential of the localized states 
for different G and compare that with the variational re- 
sults. Finally, we test the stability of the localized state 
with G = 10. While obtaining the localized state, dur- 
ing real-time propagation we suddenly change G from 10 
to 11, making the system more repulsive. Consequently, 
the root mean square size w increases and the system 
executes breathing oscillation. In Fig. [7] (d) we plot the 
root mean square size of the system during time evo- 
lution when G was increased at t = 25. The sustained 
oscillation of the system demonstrates the stability of the 
localized state. The variational equation (11) gives the 



largest value of G for which there is a localized state. 
Real solutions for the width are allowed for G < 230.7, 
beyond which there is no localization. 



tentials can be created in laboratory and such perturbed 
lattices can easily be used for BEC study. Usually a 
bichromatic lattice is used in BEC studies, where the 
periodicity of the principal OL is removed by adding a 
secondary OL of incommensurate wave length. 

Using the weakly perturbed OL we studied the prob- 
lem of localization of a noninteracting or weakly inter- 
acting BEC in ID and 3D systems using numerical and 
variational approaches to the GP equation. We studied 
the density profile and chemical potential of the localized 
states, where good agreement between variational and 
numerical results was established. The stability of the lo- 
calized state was established by the Vakhitov-Kolokolov 
criterion and complemented by numerical studies of sta- 
ble oscillation of the localized state upon small perturba- 
tion. 

In this study we were mostly concerned with station- 
ary localized states of BEC in ID and 3D in nonperiodic 
potential and their stability. We are not concerned with 
the very long-time dynamics of an arbitrary initial state 
in such a potential. For example, in a ID weakly nonlin- 
ear disordered lattice Anderson localization is destroyed 
and the field spreads sub diffusively far beyond the lin- 
ear localization zone at large times [18]. Also, in 3D if a 
wave packet propagates in a linear disordered potential, 
at intermediate times there is a diffusion regime (due to 
atoms that have energies above the mobility edge), and 
at long times, if the disorder is strong enough, a remain- 
ing localized part (due to atoms that have energies below 
the mobility edge) will emerge [4] [29]. Nevertheless, we 
demonstrate here using a mean-field model that station- 
ary localized states are possible in a disordered potential 
in ID and 3D for a relatively large nonlinearity or atomic 
interaction and this could motivate new experiments in 
the localization of BEC. 



IV. SUMMARY AND DISCUSSION 
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